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ABSTRACT 

We calculate the flux and spectrum of synchrotron radiation produced by 
high energy electrons and positrons (e^) in an advection dominated accretion 
flow (ADAF) around a black hole. The e ± are from the decay of charged 
pions which are created through proton-proton collisions. We consider both 
a thermal and power-law energy distribution of protons, and show that the 
resulting e ± synchrotron emission produces a characteristic spectrum be¬ 
tween radio and X-ray frequencies. While previous signatures of the hot 
protons were only possible at gamma-ray energies, via the production of 
gamma-rays through neutral pion decay, the present results provide a more 
observationally tractable way of probing the proton energy distribution and 
the two temperature structure in these accretion flows. We discuss a number 
of strong observational predictions from these systems, as well as the recent 
results of Mahadevan (1998) which appear to confirm the two temperature 
structure in ADAFs. We show that the results provide support for both a 
power-law and thermal distribution of protons, with at least a third of the 
viscous energy going into the power-law. 
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mechanisms: nonthermal - acceleration of particles - Galaxy: centre - radio 
continuum: galaxies 


1. Introduction 

In an advection-dominated accretion flow (ADAF), all the viscous energy generated 
is stored in the accreting gas and advected into the central star (Ichimaru 1977; Rees et 
ah 1982; Abramowicz et al. 1998; Narayan & Yi 1994,1995ab; Abramowicz et ah 1995; 
see Narayan, Mahadevan & Quataert 1998a for a recent review). If the central object 
happens to be a black hole, the advected energy is lost through the event horizon, and 
the accreting system appears very dim. Most of the viscous heating is assumed to mainly 
affect the ions, the more massive species, while the radiation is produced primarily by 
the electrons. Since the ions transfer only small fraction of their energy to the electrons 
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via Coulomb collisions, the radiative efficiency of an ADAF is much less than the total 
energy released during accretion (Rees et al. 1982). The gas in an ADAF forms a two 
temperature structure (Shapiro, Light man, & Eardley 1976; Rees et al. 1982), with the 
ions attaining nearly virial temperatures close to the black hole {T % ~ 10 12 K), and the 
electrons being much cooler ( T e ~ 10 9,5 K). 

The success of ADAF models lies in their ability to accurately predict the observed 
spectrum from a number of accreting stellar mass and supermassive black hole systems 
(see Narayan et al. 1998a for a review). The spectrum depends sensitively on the as¬ 
sumed energy distribution of the protons and electrons as well as their individual cooling 
mechanisms. The protons cool by creating neutral pions, which subsequently decay into 
gamma-rays (Mahadevan, Narayan, & Krolik 1997), while the electrons cool by various 
optically thin processes such as cyclo-synchrotron, inverse Compton and bremsstrahlung 
radiation (see eg. Narayan & Yi 1995b; Mahadevan 1997). The electrons form a thermal 
distribution (Mahadevan & Quataert 1998), and give rise to a spectrum ranging from 
radio to hard X-ray frequencies (10 9 Hz ^ v ^ 10 21 Hz), but the energy distribution 
of the protons, and therefore the resulting gamma-ray spectrum, is still unknown and 
depends on the viscous heating mechanism involved (Mahadevan & Quataert 1997; Ma¬ 
hadevan et al. 1997). If viscosity heats the protons into a thermal distribution, the 
resulting gamma-ray spectrum is sharply peaked around energies ~ 70 MeV, while a 
power-law distribution results in a power-law gamma-ray spectrum which extends to 
very high energies (Mahadevan et al. 1997). 

It is the purpose of this paper to probe the two temperature paradigm further and 
attempt to resolve what the distribution function of the protons is likely to be. We are 
therefore interested in other possible emission processes from the protons from which 
accurate spectra can be determined. In particular, we focus on the production of charged 
pions which subsequentially decay into electrons and positrons. These particles interact 
with the equipartition magnetic fields in the accreting gas, and produce synchrotron 
radiation. Although this energy leaves the plasma at lower frequencies than gamma- 
rays, it provides the exciting possibility that current sensitive and high resolution radio 
to X-ray telescopes might be used to probe the much higher energy proton distribution 
function and viscous heating mechanisms in these plasmas. 

We begin in §2 by describing all the physical processes involved in calculating the 
e ± synchrotron spectra from both a thermal and power-law parent proton distribution. 
These calculations depend only on the microphysics of particle collisions, decays, and 
emission processes.§3 introduces the physical properties of an ADAF and, using the 
results of §2, determines the general properties of the resulting e ± spectrum from an 
ADAF. §4 applies these results to solar mass and supermassive accreting black holes. 
In §5 we conclude by discussing the implications of these results on the basic assump¬ 
tions in ADAFs, and propose future observational tests which would lead to a better 
understanding of these systems. 
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2. Physical Processes 
2.1. Production of Electrons and Positrons 


The production of electrons and positrons through the decay of charged pions is 
a three step process and has been studied by various authors (see eg. Ginzburg & 
Syrovatskii 1964; Ramaty & Lingenfclter 1966ab, 1968; Stecker 1971; Bhadwar et al. 
1977; Dernier 1986ab). The colliding protons produce a charged pion, which decays 
into a charged muon, /i ± , and a muon neutrino, The charged muon decays into an 
electron, neutrino and their antiparticles. The reactions are 


p + p -> 71 ^ + A", 7r =t -> //* + 


-> e ± + v e (v e ) + (!) 

where A" represents all other decay products. The number of charged pions created per 
unit volume, per unit time, with energy E n ±, is given by (Landau & Lifshitz 1975, §12; 
cf. Mahadevan et al. 1997, eq. [4]) 


R(E n ±) = 


dN^ 


dE n ± dV 


c 

2 


, f°° , f 1 i da(E n ±; 7i,7 2 ,cos0) 

«7i / «72 / ucosft -—- 

J 7 J-i dE n ± 


x n 7l n 72 V (A - f3 2 ) 2 - (Pi x /5 2 ) 2 s 


-i 


( 2 ) 


where (/3i, /3 2 ), ( 71 , 72), (n 7i , n 72 ) are the velocity parameters, Lorentz factors, and densi¬ 
ties ofthe two colliding particles respectively, cosd = PvP 2 /\Pi\\Pi\i da(E n ±; 71 , 7 2 , cos9)/dE w ± 
is the normalized differential cross section for the production of a charged pion with 
energy E n ±. and c is the velocity of light. Following Mahadevan et al. (1997), the 
differential cross section is determined by using the isobar model for energies ^ 4 GeV 
(Lindenbaum & Sternheimer 1957; Stecker 1971; Dernier 1986a), and the scaling model 
for energies ^ 8GeV (Badhwar et al. 1977; Stephens & Bhadwar 1981; Dernier 1986b). 

At intermediate energies the differential cross section is determined by interpolating 
smoothly (cubic interpolation) between the two regimes. 

The spectrum of electrons and positrons, R(e ± ), is derived from the tt ± spectrum 
by considering the reactions, 


7T ± -> fj* + !/„(!/„), fj^ -> e ± + V e (v e ) + V^n). 

Ideally, determining the energy distribution of e ± in the laboratory frame requires three 
Lorentz transformations: one to rest frame of the 7T =1= , to determine the decay spectrum 
of the /r ± , the second to the rest frame of the /i ± , to determine the decay spectrum of 
the e ± , and the third a transformation back to the laboratory frame. However, since the 
Lorentz factor of the / u ± is small in the rest frame of the 7r ± , 

m l± + m «± 

7m± = 7^-— 

2m 7r ±m /i ± 




1.046, 


( 3 ) 
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the difference between its velocity and that of the parent 7 r ± can be neglected. In 
this approximation, R(E W ±) ~ and determining the e ± spectrum simplifies to 

requiring only two Lorentz transformations. R(E e ±) can now be obtained by performing 
the double integral (Ginzburg & Syrovatskii 1964), 


R(E e ±) = 


p' 


■R{E, 


dP' 

dEP 


dE^± dE'j 


(4) 


where the limits of integration are such that 


E e ± E e ± | ^ p e ± r Yi T ±(d n ±. 


Here, dP'/dE' e ± is the probability of creating an e ± with energy E' e ± and momentum 
p' e ± in the rest frame of the decaying muon, and is given by (Ginzburg & Syrovatskii 
1964), 

dP' _ 16 E' e± , 
dE' e± ~ ml Pe 


3-4^i 


m„ 


(5) 


2.2. Synchrotron Emission 


If the plasma in which the e^’s are created has a magnetic held, the particles inter¬ 
act with the magnetic held and cool by radiating synchrotron emission. The synchrotron 
radiation produces a characteristic spectrum which depends directly on the energy dis¬ 
tribution of the e . Since the distribution of e ± is determined by the proton energy 
distribution (cf. eqs. [||, H), the resulting synchrotron spectrum provides a probe 
to the energy distribution of the protons. Here, we provide the basic equations that 
determine the e ± synchrotron spectra, and show (§2.3) how different initial proton dis¬ 
tributions give rise to different spectra. In what follows all quantities refer to the created 
e 1 * 1 , unless otherwise noted. 


Calculating the synchrotron emissivity, L s v {v), requires evaluating 

r'J max 

L s u du= / IV( 7 ) j(z/, 7 ) d'y erg s ” 1 cnr 3 Hz -1 , ( 6 ) 

J 'Tmin 

where 7 is the Lorentz factor of the created e ± 1 j(u, 7 ) is the synchrotron emission 
formula averaged over all particle directions, and A^(y) is the steady state distribution 
of e ± . Since the e ± created are highly relativistic (§2.3, Fig. |1]), and extend to very high 
energies, we can set 7 min ^ 1, and 7 max —> 00 . In this limit, j(z/, 7 ) takes the relativistic 
form (eg. Rybicki & Lightman 1979), 

jO, 7) = So F , 


So = 


47TV^r 2 


e u B 


eB 


— 2 T Vb i 


u b = 


3c 


v. 


2irm e c 


(7) 
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where we have averaged over all particle directions. Here F(x) is an integral over mod¬ 
ified Bessel functions (Rybicki & Lightman 1979, eq.[6.31c]), S 0 = 1.56 x 10 -22 B, and 
vb = 2.8 x 10 6 Hz is the cyclotron frequency. Integrating equation (j7|) over all frequencies 
gives the total energy radiated per particle, 

poo 

Es= j(v, 7 ) dv = 1.06 x 1CT 15 y 2 B 2 ergs -1 . ( 8 ) 

Jo 

For large 7 , the e ± radiate their energy very efficiently, and their steady state distri¬ 
bution, N( 7 ), is determined by the competing effects of the “creation” and “depletion” 
of particles. At a given energy E, the colliding protons produce R(E) electrons and 
positrons. These e, however, cool very efficiently, thereby reducing the number of e ± 
at energy E. The steady state distribution can therefore be determined by solving the 
Fokker-Planck equation without the stochastic term, 

^fi 2 l = -A | 7i v( 7,« ) ]+fl ( 7,i ) . P) 

with dN( 7 , t)/dt = 0. Here, 7 = Es/m e c 2 is the cooling rate of the and R( 7 ) is their 
injection rate (cf. eq. [Q). I 11 steady state, the equation requires the net flux of particles 
between 7 x and 72 to be equal to the rate of injection of particles between these two 
Lorentz factors. Alternatively, the equation can be thought of as a particle conservation 
equation in energy space with 7 being the velocity along the energy axis. 

Setting dN{^,t)/dt = 0, equation (g) gives 

iV(7) = -ic(7), 

7 

poo 

c( 7 )=/ 

J'l 

where C( 7 ) is the total number of injected e ± above Lorentz factor 7 
(|6j), the resulting synchrotron spectrum can be written as 

L s u dv= f C( 7 ) ^ V ] ^ d'y erg cm -3 s -1 Hz -1 . (11) 

•Amin 7 


( 10 ) 

Using equation 


2.3. Synchrotron Spectra from Electrons and Positrons. 

Using equations (§, ©, and flll]) , the synchrotron spectra from the e ± can be 
determined, and depends only on the initial proton energy distribution. We consider here 
two extreme distributions: a relativistic Maxwcll-Boltzmann, N(E) oc E 2 exp(— E/kT), 
and a power-law with energy index s, N(E) oc E~ s . In almost all areas of astrophysics, 
the true proton distribution is probably some combination of the two. 
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Before determining the exact synchrotron spectra, we point out two general prop¬ 
erties of equation (|Tl]). First, the total synchrotron power, P total , is given by 


-Ptotal 


L s v dv = m e c 2 / ( 7 ( 7 ) dy erg s 


-1 


( 12 ) 


which is equal to the total energy of e ± that are created per second. This depends 
only on the injected energy spectrum and is independent of the synchrotron emissivity, 
j(u, 7 ), and therefore of the magnetic field. We might expect that since particles in a 
very low magnetic field radiate synchrotron emission inefficiently (cf. eq.||), the total 
synchrotron emission from this plasma would be low. ffowever, since we are interested in 
steady state distributions, the decrease in the amount of cooling is exactly compensated 
by an increase in the steady state number density (cf. eq. »■ Similarly, higher 


magnetic fields lead to more efficient cooling, which results in a lower steady state 
number density of particles. The two effects exactly compensate for each other and 
leaves the total cooling rate unchanged. 


Second, the synchrotron spectrum, in a uL u versus v plot, is only a function of the 
dimensionless frequency x = v / v b (cf- eqs. 10)- This shows that while the intrinsic 
shape of the synchrotron spectrum does not change as the magnetic field changes, the 
position of the spectrum does move horizontally to the right (left) as the magnetic field 
is increased (decreased). 


2.3.1. Thermal Proton Distribution 

Using a normalized relativistic Maxwell-Boltzmann distribution for n 71 , n 72 in equa¬ 
tion (j2[), RtyEifi) is obtained by performing the integrals given in Dernier (1986a; see also 
Mahadevan et al. 1997, eq. [12] with the appropriate changes for 7r ± cross-sections). The 
results are then used in equation ([|) to determine the resulting production spectrum, 
R(E e ±), of electrons and positrons. 

Figure |l]a shows R(E e ±) for different values of the dimensionless proton temperature 
9 P = kT p /m p c 2 = 0.05,0.1,0.2. The spectra rise at low energies, peak near E e ± ~ 
33 MeV, and then decline exponentially at higher energies. 

An important feature shown in Figure |I]a is the effect of small changes in tem¬ 
perature on R(E e ±). At low temperatures, only protons in the exponential tail of the 
Maxwell-Boltzmann distribution have energies above threshold to produce pions and 
therefore electrons and positrons. Increasing the temperature by a factor of two from 
0.05 to 0.1, dramatically increases the number of protons with energies above threshold, 
thereby increasing R(E e ±) by more than an order of magnitude. Changing the temper¬ 
ature again, from 0.1 to 0.2, has a less drastic effect on the production rate. At these 
high temperatures the protons responsible for most of the 7r ± production are no longer 
in the tail of the Maxwell-Boltzmann distribution, but rather have energies of order the 
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average energy of the plasma. Therefore changing the temperature by a factor of two, 
has the effect of increasing R(E e ±) by nearly the same amount. 

Figure |l]b shows the resulting steady state distribution of (cf. eq. [p~G|1 ). At 
low energies, the product 7 ^( 7 ) is a constant, which gives a power-law distribution 
TV( 7 ) oc 7 -2 , while at higher energies N( 7 ) decreases exponentially. The low energy 
power-law is the result of high energy e ± that lose their energy due to efficient synchroton 
cooling. 

Figure |l]c shows the resulting e ± synchrotron spectrum for three values of the mag¬ 
netic field B = 10, 10 4 , 10 8 Gauss (solid, dashed, and dotted lines respectively), and for 
6 P = 0.2 (lower values of 9 P result in identical spectra, but with much lower fluxes). As 
expected, increasing the magnetic field merely shifts the spectrum to higher frequencies, 
and leaves the total flux and shape of the spectrum unchanged. The change in the steady 
state distribution A r (y) from a power-law at low energies (oc 7 -2 ) to an exponential at 
high energies, results in a synchrotron spectrum that rises as vL v oc i/ +0 ' 5 , and then 
turns over exponentially (Rybicki & Lightman 1979). The frequency, u t , at which the 
spectrum turns over is determined by the Lorentz factor of the e ± responsible for most 
of the radiation (7 ~ 100, see Fig. 0 b), as well as the magnetic field strength, and is 
given by v t ~ 3 x 10 6 y 2 B ~ 10 10 B. 

2.3.2. Power-law Proton Distribution 

Using a normalized power-law proton distribution, n 7 = (s — l)y _s , with energy 
index s, equations ([?[) and (|j) are used to determine the production rate, R(E e ±), of 
electrons and positrons. The calculation is identical to that of 7T° production with the 
appropriate changes for 7 r ± cross-sections (see Mahadevan et al. 1997). 

Figure |l]d shows R(E e ±) for three values of the energy index s. The spectrum rises 
at low energies, turns over at E e ± ~ 35 MeV, and extends as a power-law, E~ s , with 
the same energy dependence as the parent proton distribution (Ginzburg & Syrovatskii 
1964). While the energy at which the spectra turn over are similar to the thermal 
case (cf. Fig. |Ija), the non-thermal spectra are much broader and extend to very high 
energies. In addition, the production rate is not very sensitive to the energy index s. An 
increase in s results in a small decrease in the number of protons above the threshold 
energy, which moderately decreases the total production rate. 

Figure [l]e shows the resulting steady state e ± distribution. Similar to the thermal 
case, the product 7 ^( 7 ) is constant at low energies which gives a power-law distribution 
W(y) oc y -2 . At higher energies, however, the shape of the distribution deviates sub¬ 
stantially from the thermal case, and extends to very high energies as A^(y) oc 7 W -1 ) 
(cf. eq. 0 ). 


The resulting synchrotron spectrum is shown in Figure |ljf for three values of the 
magnetic field B = 10, 10 4 , 10 8 Gauss, and for s = 2.75. As in the thermal case, increas¬ 
ing the magnetic field merely shifts the spectrum to higher frequencies and leaves the 
total flux and shape of the spectrum unchanged. The change in the index of the steady 
state distribution from low to high energies, results in a synchrotron spectrum that rises 
as uL v oc z/ +0 - 5 , which then steepens to uL v oc z/ _S//2 (Rybicki & Lightman 1979). There¬ 
fore, while different values for the proton energy index s result in identical spectra at low 
frequencies, smaller values of s give rise to harder spectra at high frequencies. For all s, 
the frequency at which the spectrum turns over is given by, ^~ 3 x 10 6 7 2 B ~ 10 11 B, 
where we have set 7 ~ 200 (cf. Fig. p). 


3. Application to ADAFs 

The results from the previous section can now be used to determine the total e ± 
synchrotron spectrum from an ADAF. Since an ADAF is well approximated by a series of 
concentric spherical shells (Narayan & Yi 1995ab), with the properties of the gas varying 
with radius, the total e ± spectrum is obtained by determining the synchrotron spectra 
from each shell, and then propagating the resulting spectra through the accretion flow 
to the observer (§3.2). Determining the spectrum from each shell requires a knowledge 
of the proton densities and magnetic field strengths at each radius. These quantities are 
obtained by solving for the global structure of an ADAF. 


3.1. ADAF Equations and Properties 

The global structure of an ADAF is characterized by four parameters: the viscosity 
parameter a (Shakura & Sunyaev 1973), the ratio of gas to total pressure /3 a d v , the mass 
of the central black hole M, and the accretion rate M. Given these parameters, the 
global structure of an ADAF can be calculated ( Gammie & Popham 1998; Popharn 
& Gammie 1998; Narayan, Kato, & Honma 1997; Chen, Abramowicz, & Lasota 1997), 
and for the results presented here, we use the new fully relativistic, self-consistent, 
global solutions in the Kerr metric developed by Popham & Gammie (1998). In general, 
ADAFs appear to have large values of a ^ 0.1, and the uncertainty in the value is a 
factor of a few (Narayan 1996). Most ADAF models in the literature use a = 0.3 (see 
Narayan et al. 1998a). For the magnetic field strength, we consider two extreme values 
of /3 a dv = (0.5, 0.9), which corresponds, respectively, to an equipartition field, and to one 
where the field contributes negligibly (~ 1/10) to the total pressure. 

Assuming that the mass fraction of hydrogen is X = 0.75, the total numberdensity 
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of protons is obtained from the continuity equation, [] and is given by, 


n v (r) = — P_J_ _ g x 1 r y ir 2 

' m p 

= m~ 1 rhh(r) cm _3 , 


v(r) 


(l- [v(r)/c} 2 \ 

c 

cos 6 * h 

{ l-^” 1 ) 


( 13 ) 


where p(r ) is the mass density and v(r ) is the radial velocity of the gas. Here, m = M /M 0 
is the mass of the black hole in solar mass units, m = M/M Edd is the accretion rate in 
Eddington units with with M Edd = 10L Edd /c 2 = 1.38 x 10 18 m g s -1 , and r = R/R s is 
the radius in Schwarzschild units where R s = 2 GM/c 2 = 2.95 x 10 5 m cm. This equation 
differs from that in Mahadevan et al. (1997) by the two last terms on the right. The 
factor (cos^h)^ 1 accounts for the slight non-spherical geometry of the flow (Narayan, 
Barret & McClintock 1997a, Appendix A), and the term in parentheses is a general 
relativistic correction factor (Popham & Gammie 1998). We obtain v{r)/c and cos$h 
from the global ADAF solutions of Popham & Gammie (1998). 


The magnetic field strength is determined by the parameter f3 aa v , and is defined by 
(Narayan & Yi 1995b), 

B 2 (r) (l~P a dv\ 

Pg > 


24 7T 


Pg = (3advp(r) c 2 s (r ) = 


/3a,dv 

p(r)k T p p{r)kT t 


(14) 


fii Trip fi e Trip 

where p g is the gas pressure, c s (r) is the isothermal sound speed, and pi = 1.23, p e = 1.14 
are the effective molecular weights of the ions and electrons respectively. Since T p T e 
we neglect the second term to obtain 

2 


T p (r ) = 1.34 x 10 13 /3 adv 


c s (r) 


K. 


(15) 


We obtain c s {r)/c from the global ADAF solutions (Popham & Gammie 1998). 


Given n p (r) and B{r) from equations (pA[) , O. and (HD, the spectrum can now 
be determined. 


3.2. ADAF Spectra 

In addition to the numberdensity of protons and magnetic held strength, the spec¬ 
trum from each radius, and therefore the total e^ synchrotron spectrum from the ADAF, 


lr The continuity equation is given by, 


( 1-Rs/R 


1/2 


M = 4irR 2 p(R)v(R) cos^h 


1 - [v(R)/c\ 2 
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also depends on the energy distribution of the protons (§2.3). Since our knowledge of 
the viscous heating of protons is unknown, it is still uncertain as to whether the process 
leads to a thermal or power-law distribution in proton energies, or possibly a combina¬ 
tion of the two. Below, we explore both distributions and calculate the e ± synchrotron 
spectrum associated with them. 

At each radius, we write the numberdensity of protons as a function of energy, 
n p (r, 7 ), as the product of a normalized velocity distribution, which depends on the 
temperature at each radius, and the total numberdensity of protons (see eg. Mahadevan 
et al. 1997), 

M r > 7 ) = n p (r)n^[y,e p (r)}, (16) 

with 

/ oo ro o x 

n 1 ['y,9 P (r)]d'y = l, ^ (7 - 1) n 7 [y, 0 p (r)] dy = - d p (r). (17) 

The total numberdensity of protons, n p (r), is determined from equation (|l3|) , and 
for the present results we consider a normalized velocity distribution, n 7 [y, 0 p (r)] which is 
either a relativistic Maxwell-Boltzmann or power-law distribution. With this functional 
form of n p (r, 7 ), equations (Q) and ([3|) can be evaluated to give the total production of 
e^, and therefore the total synchrotron emission. 

Determining the synchrotron spectrum from the e ± in an ADAF, requires three 
additional pieces of physics to those presented in §2. First, the emergent spectrum at 
radii close to the black hole must be corrected for gravitational redshift effects. This is 
included in the calculations that follow, and has the effect that L s v (y) observed at large 
radii is now L s v [v{l — 1/r 1 / 2 )]. This shifts L s v redward by the gravitational redshift factor. 
The inclusion of gravitational redshift will change the entire optically thin synchrotron 
spectrum determined in §2.3. However, this will only affect the spectrum dramatically 
within a few Schwarzschild radii. 

Second, in addition to calculating the optically thin synchrotron spectrum at each 
radius in the flow (§2.3), absorption processes also need to be included. This will change 
the resulting synchrotron spectrum at frequencies where the plasma is optically thick. 
In particular, since the thermal electrons in an ADAF radiate highly self-absorbed syn¬ 
chrotron emission (Narayan & Yi 1995b; see also Mahadevan 1997), these electrons will 
absorb any radiation that is above the local black body value. Therefore, to account for 
this absorption process we replace any “excess” synchrotron emission from the e ± with 
the local black body spectrum at each radius. In a highly self-absorbed plasma, this is 
a good approximation to the exact radiative transfer calculations. 

In addition to the absorption process described above, self-absorption by the created 
e ± is also included. However, in the calculations that follow, we find that absorption 
from the thermal electrons is always more important. Since the amount of absorption 
depends on the physical size of the emitting region, the calculations for the rest of this 
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section are for a black hole mass m = 10 6 , accreting at a rate m = 10 4 . §4 explores the 
results for different black hole masses and accretion rates. 


Finally, the calculations and spectra discussed in §2 assumes that a steady state dis¬ 
tribution has been established. The only timescale of interest there was the synchrotron 
cooling time, ffowever, in an ADAF, the accretion time, t ac , is another important 
timescale in the problem. If the synchrotron cooling time is longer than the accretion 
time, a steady state distribution will not be able to form, and while the e ± will still 
radiate synchrotron emission, the spectral results of §2 are no longer valid. Setting 
the synchrotron cooling time t s = 7 /% — 7.4 x 10 8 / 7 H 2 , equal to the accretion time 
t ac — 1-8 x !CT 5 cY 1 mr 3 / 2 (Mahadevan & Quataert 1998), gives 


m ^ 4 x 10 4 




/ 1 ~ Aidv \ 


(18) 


where we have used the self-similar scaling of the magnetic field with radius (Narayan 
& Yi 1995b; see also Mahadevan 1997), and have set 7 ~ 200, the Lorentz factor of the 
e ± responsible for most of the radiation (cf. Fig. [l]e) . For systems where m ^ 10~ 4 , 
the timescale for the e ± to cool via synchrotron radiation is less than the accretion time, 
and a steady state distribution will be established. The results of §2 are therefore valid. 
For lower accretion rates, the spectrum will be modified, but we do not consider this 
regime since the luminosities from these systems will be too low to be of observational 
consequence. All published models of ADAFs in the literature have m > ICY 4 . 


3.2.1. Thermal Distribution 


Using equation (p~3|) the synchrotron emissivity per unit scaled volume in an ADAF 
is given by 

dE 


L v = 


= 2.57 x 10 


16 mm 2 h 2 (r) 


LI erg s Hz 


,-1 


(19) 


(cf. eq. |ll|), and dV = 


dt dV dv 

where L s u is the synchrotron emissivity for n p = 1 cm -3 
(2.95 x 10 5 ) 3 m 3 dV. Here, dV = 47 rr 2 dr, is the scaled volume in Schwarzschild units. 
The dependence of L v on the energy distribution of the protons is taken into account by 
the form of R(E e ±) that is used to evaluate L s u (cf. Figure [L]a). As expected, equation 
(P~9|) shows that the luminosity of the optically thin part of the spectrum increases as 
mrh 2 similar to the 7 -ray spectrum (Mahadevan et al. 1997). However, this scaling is 
no longer valid for frequencies where the spectrum is self-absorbed. 


All the emission is produced in the inner most regions of the flow (r ^ 5) where the 
proton temperatures are high enough to produce pions. At larger radii, the temperatures 
drop rapidly, and almost none of the protons have energies above the threshold for 
pion production (cf. § 2 ). Since the synchrotron emission originates from a few 
Schwarzschild radii, all the self-absorption is from the thermal electrons at r ^ 5. 
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Figure 2a shows the results for two values of the viscosity parameter a (0.1 and 
0.3) and two values of /? adv (0.5 and 0.9). At relatively high frequencies, the spectra are 
optically thin and show the characteristic rise ( vL v oc u +0 - 5 ) and exponential turnover 
expected from a thermal proton distribution (§2, Fig. 0 c). At lower frequencies, however, 
the emission is highly self-absorbed by the ambient thermal electrons. The resultant self- 
absorbed spectrum is Rayleigh-Jeans and has a spectral dependence uL u oc zA Although 
the distinct spectral change offers an interesting testable prediction, the luminosity from 
this process is negligible compared with the synchrotron luminosity from the thermal 
electrons at these frequencies. The thermal e ± spectrum is therefore undetectable (see 

§ 4 ). 


For a given a, increasing f3 adv leads to an increase in gas pressure, which therefore 
increases the proton temperature (cf. eq.JT5[). Changing /9 adv from 0.5 to 0.9, causes 
the temperature to increase by a factor ~ 2 , and since the flux is extremely sensitive 
to temperature (cf. § 2 ), the synchrotron flux increases by nearly 2 orders of magnitude. 
For a fixed f3 adv , changing a also affects the synchrotron luminosity. Since tt ± production 
scales as n 2 oc a~ 2 (Narayan & Yi 1995b), we might expect that decreasing a increases 
the total luminosity. In the case of high /3 adv this is the case, but in the more detailed 
global solutions in the Kerr metric, the slight change in temperature with a in the in¬ 
nermost regions of the flow, for low /3 adv , compensates for the increase in numberdensity, 
and the net result is a decrease in the total flux. 


3.2.2. Power-Law Distribution 

For a power-law distribution of protons, the total number of protons is still given 
by equation (pA|) , but the distribution n 7 [ 7 , 6 * p (r)] is no longer Maxwellian. Following 
Mahadevan et al. (1997) we write the energy distribution of protons as 

nb,9 P (r)\d'y= {[1 - COO] 5(7 - 1) + (s - 1) COO 7' s } <*7, ( 20 ) 

where a fraction 1—CO") of the protons have 7 ~ 1 , and a fraction COO are in a power-law 
tail with index s. The fraction C( r ) is fixed by the energy requirement (cf. eq. ||17|| ): 

CM = \ (s - 2 ) ( 21 ) 

Using equation (|20|) in equation (Q) gives the total charged pion spectrum. In an ADAF, 
(•Cl, and Mahadevan et al. (1997) have shown that most of the contribution to the 
pion spectrum comes from protons in the power-law tail colliding with protons at rest. 
Equation (^) can therefore be simplified to 

n [7>0pOO] d 7- 0 -lRM 7“ s ^7, 

and the e ± synchrotron emissivity per unit scaled volume takes the form 

Lu = ofnDn = 2 ' 57 x 1Ql6 mrh2 fi ?( r ) C O’) t 1 ~ C O’)] K, ( 22 ) 

at a V du 
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where L s v represents the synchrotron emissivity using a proton distribution n( 7 ) = (s — 
l) 7 _s (cf. eq.jnj] and Fig. |l]f). 


Figure 2b shows the resulting spectrum for a proton energy distribution with power- 
law index s = 2.75 for different values of a = (0.1, 0.3) and /d a d v = (0.5, 0.9). Similar to 
the thermal case, the optically thin part of the spectrum scales as mth 2 (cf. eq. |[ 22 || ), but 
at lower frequencies this scaling is no longer valid since the spectrum is self-absorbed. 
For a fixed a increasing /3 a d v decreases the magnetic field which shifts the spectrum to 
lower frequencies (§2.3, Fig. Iff). For fixed /3 a d v decreasing a increases the numberdensity 
(n v oc a -1 ) which increases the optically thin flux. 


Unlike the thermal case, the emission comes from nearly all radii in the ADAF. In 
particular, since the protons are assumed to be accelerated into a power-law at all radii, 
most of the high frequency emission originates from r ^ 30 while the low frequency 
emission comes from r ^ 10 4 . To understand this effect, recall that an e^ with Lorentz 
factor 7 radiates most of its energy at frequencies v oc q 2 B. Since 7 is fixed at all radii 
by the microphysics of particle decays (§2 and Fig. Iff), the frequency at which most of 
the synchrotron radiation emerges depends only on the magnetic field. As the magnetic 
field decreases with increasing radius, lower frequency emission occurs farther away from 
the black hole. 


At relatively high frequencies, the spectra are optically thin and show the charac¬ 
teristic rise (vL v oc z/ +a5 ), turnover, and power-law tail {uL v oc u~ s ^ 2 ) expected from a 
power-law proton distribution (§2, Fig. Iff). At lower frequencies, however, the emis¬ 
sion is highly self-absorbed by the ambient thermal electrons, and predicted spectrum 
is determined by self-absorption at all radii (§3.2). We show below (§§4 and 5) that the 
expected spectral shape provides an interesting prediction for the ADAF paradigm, and 
in the case of our Galactic Centre, Mahadevan (1998) has shown that these predictions 
do indeed agree quite well with the observations. 


4. Results: Application to Solar and Supermassive Black Holes 

This section combines the spectra produced by both the protons and electrons in 
an ADAF. The bulk of the electrons cool via thermal synchrotron, bremsstrahlung and 
inverse Comptonization (Narayan & Yi 1995b; see also Mahadevan 1997), while the 
protons cool via the production of neutral and charged pions. In the results that follow, 
the 7 -ray spectrum, produced by the decay of neutral pions, has been updated from 
Mahadevan et al. (1997) to include general relativistic effects, gravitational redshift, 
and the non-spherical geometry of the flow (§3). 
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4.1. Thermal Distribution 

For a thermal distribution of protons, Figure |3]a shows the total spectrum from 
ADAFs around various black hole systems corresponding to m = (10,10 6 ,10 9 ) and 
m = (10 - 2 ,10 -4 ). The solid, dashed and dotted lines correspond respectively to the 
total emission from the ADAF, the emission from the bulk of the thermal electrons, and 
the emission from the eV As expected, the low proton temperatures in an ADAF do 
not allow for a significant amount of e ± synchrotron emission, and the luminosities are 
too low to be of observational interest. 


4.2. Power—Law Distribution 

For a power-law distribution of protons, Figure |3p:i shows the total spectrum from 
ADAFs around the same black hole systems corresponding to m = (10,10 6 ,10 9 ) and 
fn = (10~ 2 ,10~ 4 ). The solid, dashed and dotted lines correspond respectively to the 
total emission from the ADAF, the emission from the bulk of the thermal electrons, and 
the emission from the e 1 * 1 , for the three different values of the proton energy index s. 
The heavy solid line corresponds to the total ADAF spectrum for s = 2.75. Unlike the 
thermal case, the power-law e ± do produce significant amounts of detectable emission. 

For all black hole masses, Figure [|b shows that bremsstrahlung emission from the 
thermal electrons has a cutoff at hv ~ kT e , which corresponds to v r^j 10 2 ° Hz, and 
inverse Compton scattering of the electrons also turns over at these frequencies due 
to the decrease in the Klcin-Nishina cross-section. The emission for v ^ 10 2 ° Hz is 
therefore only due to proton cooling. All spectra show that the e ± naturally provide 
a hard power-law tail which extends up to ~ 100 MeV. Beyond these energies, the 
decay of neutral pions into 7 -rays dominate the spectrum, and produce a distinct rise 
in the spectrum which turns over at a few GeV and extends as a power-law to very high 
energies. An interesting consequence is that although the spectral index of the MeV 
and GeV spectrum are different, they are not independent of each other. For a parent 
proton distribution with energy index s, the MeV and GeV spectrum have spectral 
dependencies vL v oc v ~ s / 2 and reS p ec tively (Mahadevan et al. 1997). 

For frequencies v ^ 10 2 ° Hz, the dominance of the proton cooling spectrum, com¬ 
pared with the emission from the thermal electrons, depends sensitively on the accretion 
rate m, and to a lesser degree on the mass m of the black hole. For low mass black holes 
(eg. Cyg X-l, A0620) the proton cooling spectrum is almost always unobservable for 
all accretion rates. For higher mass black holes, however, the proton cooling emission is 
observable for v ^ 10 11 Hz and between 10 14 — 10 1 ' Hz, and depends crucially on m. 

For supermassive black holes, Figure |3|b shows that between frequencies 10 14 — 10 1 ' 
Hz, the turn over in the e ± synchrotron spectrum, is observable only for low accretion rate 
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systems. Increasing the accretion rate, increases the amount of thermal inverse Compton 
scattering (oc rh 3 ; Narayan & Yi 1995b; Mahadevan 1997), which increases faster than 
the amount of synchrotron radiation (oc m 2 ). The luminosity from thermal inverse 
Compton scattering will therefore dominate the spectrum at these frequencies. The 
e ± spectrum is therefore unobservable for high accretion rate AGN such as NGC 4258 
(m ~ 0.01; Lasota et al. 1996), while the spectrum dominates the infrared to soft X-ray 
band for low accretion rate systems like our Galactic centre (rh ~ 10“ 4 , Narayan et al. 
1998b; Mahadevan 1998). 

For frequencies v A 10 11 Hz, the radio spectrum from ADAFs around supermassive 
black holes is substantially modified by the radiating e ± . Figure |3|b shows that the 
additional synchrotron emission has two effects: ( 1 ) it increases the radio flux by as 
much as an order of magnitude at certain frequencies, and ( 2 ) changes the shape of the 
spectrum by producing an observable spectral break. As discussed in §5, Mahadevan 
(1998) has shown that the observed spectral break from radio observations of the Galactic 
centre (Sgr A*) can be explained quite well by e ± synchrotron radiation. 

It is interesting that the synchrotron radiation from the e ± actually produce more 
emission at low radio frequencies than the thermal electrons. Since the thermal electrons 
radiate highly self-absorbed synchrotron emission, the emission at each frequency in the 
radio, Uf, corresponds to black body emission at a definite radius, rf (Mahadevan 1997). 
Any “excess” emission that is greater than the local black body value will be highly 
self-absorbed. In this case, how does the inclusion of synchrotron radiation from the 
e ± increase the radio emission? The apparent paradox is resolved by noting that the 
enhanced emission at Vf results from the high energy e ± radiating synchrotron emission 
at larger radii, r > r/. While any excess radiation emitted at r < 77 is completely 
self-absorbed by the thermal electrons at 77 , the e 1 * 1 at radii r > 77 radiate most of their 
energy at v ~ 7 2 z/g ~ 17 , which is now optically thin to the local black body, and free 
to escape from the plasma. The excess radio emission therefore comes from high energy 
e ± radiating at larger radii. 


5. Discussion and Conclusions 

ADAFs are the only hot accretion flow models that attempt to self-consistently 
incorporate the viscous hydrodynamics, thermal structure, and radiative processes in the 
inflowing gas. The models have been worked out in considerable detail, and the success 
of ADAFs therefore lies in their ability to make robust observational predictions. The 
numerous applications of ADAFs to solar and supermassive black holes (eg. Di Matteo 
et al. 1998; Esin et al. 1997, 1998; Fabian & Rees 1995; Reynolds et al. 1996; Narayan 
et al. 1995; see Narayan et al. 1998a for a review) have not only led us to a better 
understanding of these systems, but have also reinforced our belief that ADAFs are a 
substantial step in the right direction, towards a possibly more detailed solution of low 
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luminosity accreting systems. 

In order for the ADAF solutions to exist two basic assumptions in plasma physics 
must be satisfied: (1) the existence of a hot two temperature plasma, and (2) the viscous 
energy generated primarily heats the protons. If either of these assumptions fail, then 
the ADAF solution, and the resulting spectra that they produce, cease to exist. To 
understand this, consider first the existence of a two temperature plasma. In an ADAF, 
all the viscous energy is stored in the protons of the gas and only a small fraction of this 
energy is transferred to the electrons to be radiated away. This fraction is determined 
by Coulomb collisions and the large temperature difference of the protons and electrons. 
If some other mechanism (eg. a plasma instability) is more efficient in coupling the 
protons and electrons to equilibrate their temperatures, then all the stored energy in 
the protons, would be transferred to the electrons to be radiated away. No energy will 
be advected with the accreting gas, the accretion flow becomes very luminous, and the 
ADAF solution does not exist. 

Second, the protons become hot and advect the viscous energy because viscosity 
is assumed to primarily heat the protons. If the viscously dissipated energy went into 
heating the electrons instead, then this energy would be immediately radiated away, 
resulting in a luminous system. Again, no energy would be advected by the gas, and the 
ADAF solution does not exist. The validity of these assumptions are therefore crucial 
for the existence of ADAFs. 

This paper has focussed on providing additional observational tests for these two 
assumptions, by presenting detailed spectra from both the protons and electrons. We 
have concentrated on the cooling of protons via the production of charged pions, which 
subsequently decay into positrons and electrons. These particles interact with the local 
magnetic fields to produce synchrotron radiation. Until recently, observational signa¬ 
tures of hot protons was relegated to very high energies, where neutral pions decay to 
produce gamma-rays (~ few GeV). The present results not only probe the two temper¬ 
ature paradigm further, but also allow us to predict the existence of hot protons in a 
more observationally tractable regime of the electromagnetic spectrum. 

It is important to understand the connection between the spectra from the protons 
and electrons in ADAFs. Once the spectrum from the thermal electrons is determined, all 
the parameters in the ADAF are fixed. Properties of the gas such as the numberdensity, 
temperatures, and magnetic held strengths, cannot be changed, and the resultant spectra 
from proton cooling depends only on the physical processes in the plasma. In particular, 
the energies of the created e ± 1 their steady state numberdensity distributions, and shape 
of the resulting synchrotron spectra (cf. Fig. [I]), depend only on the microphysics of 
particle collisions, decays, and radiative processes. It is this strong interplay between 
the proton and electron spectra that allows testable predictions of the basic assumptions 
of ADAFs, and also allows a probe into the proton energy distribution in these hows. 
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Since it is not understood whether viscous heating produces a thermal or power- 
law distribution of proton energies, we have calculated the e ± spectrum due to both 
distributions; the true distribution will be some combination of the two, and can be 
determined by taking the weighted sum of the thermal and power-law spectra. Recently, 
there have been a number of theoretical investigations which consider particle heating in 
ADAFs (Bisnovatyi-Kogan & Lovelace 1997; Blackman 1998; Quataert 1998; Gruzinov 
1998; Quataert & Gruzinov 1998; see Narayan et al. 1998a for a review), some of which 
also argue for and against both proton distributions. Here, we have attempted to answer 
this question observationally by providing accurate spectra from these distributions. 

If the proton distribution is thermal, the resulting synchrotron spectrum peaks 
between frequencies 10 12 — 10 1 ' Hz, and produces very low luminosities. In addition, 
the radiation by the thermal electrons dominates the emission from the radio to hard 
X-rays, which results in an e ± spectrum that is completely unobservable (Fig. |3ja). In 
this case, the only observable signature of hot protons would be in the gamma-rays, 
but the predicted fluxes are much too low to be detected even by the next generation 
gamma-ray telescope, GLAST. 

However, if the protons have a power-law distribution of energies, the e ± syn¬ 
chrotron spectrum can be observed. For supermassive black holes, the emission appears 
in the radio, and, depending on the accretion rate, also between infrared and soft X-ray 
energies. In particular, the best systems to look for such radiation would be in super- 
massive black holes with low accretion rates. Fortunately, the supermassive black hole at 
our Galactic centre (Sgr A*) offers such an opportunity. It is the closest low luminosity 
AGN that we know, and has been observed quite extensively. 

Recently, Narayan et al. (1998b) have used an ADAF model to explain the spectrum 
and low luminosity from Sgr A*. For the canonical values of a = 0.3 and /3 = 0.5, and 
the dynamically measured mass M = 2.5 x 10 6 M 0 (Haller et al. 1996; Eckart & Genzel 
1997), they vary just one free parameter, rh, to fit the X-ray flux, and show that the 
resulting spectrum from the thermal electrons explains the low luminosity and observed 
broad band spectrum quite well. The model, however, has difficulties in explaining the 
non-uniform radio spectrum. It is important to realize that once the thermal electron 
spectrum is fixed, all the parameters in an ADAF are determined, and there are no more 
free parameters available. The non-uniform radio spectrum has therefore always been 
unexplained by the ADAF model. 

This discrepancy has recently been solved by considering e ± synchrotron emission 
from the ADAF around Sgr A*(Mahadevan 1998). Without changing any parameters, 
and just including another physical process - synchrotron emission from a power-law 
proton distribution, Mahadevan (1998) shows that the ADAF model naturally repro¬ 
duces the observed radio spectral break at ~ 86 GHz, and accounts for the low frequency 
luminosities. This is shown by the solid line in Figure which provides quite compelling 
evidence that a two temperature plasma probably does exist in an ADAF around Sgr 
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A*. More importantly, the results provide, for the first time, observational support for 
the basic assumptions in an ADAF (Mahadevan 1998). 

While one of the ADAF assumptions requires almost all of the viscous energy to 
be transferred to the protons, it does not indicate whether this energy should heat the 
protons into a power-law or thermal distribution. This depends on the details of viscous 
heating, and it is interesting that the recent results from Sgr A*allows an observational 
answer to this question. If we characterize the amount of viscous energy that goes into 
a power-law by A, then the case when no energy goes into the power-law distribution 
is A = 0, which is shown in Figure [|a; A = 1 corresponds to Figure |3|b. The question 
arises as to what the most likely values of A can be, ie. how much of the viscous heating 
is transferred to a power-law distribution and how much to a thermal one. To assess 
this, we use the results of Mahadevan (1998) as our baseline model. A = 1 corresponds 
to the solid line in Figure |3j. while A = 0 corresponds to the long dashed line (Narayan 
et al. 1998b). The dotted, short dashed, and dot dashed lines correspond to spectra for 
A = 0.5,0.25, and 0.1 respectively. 

Figure fj] shows that the results are not very sensitive to A; this factor will only 
appear linearly in equation ([A]). While the luminosity of the optically thin part of the 
power-law spectrum decreases linearly with A, the self-absorbed radio part is less sensi¬ 
tive. For A ^ 0.3, the agreement with the radio spectrum does not change significantly. 
This implies that the observations are consistent with approximately a third or more 
of the viscous energy going into a power-law distribution of protons, and the rest into 
a thermal or quasi-thermal one. The results therefore provide clues to the physics of 
viscous heating in these flows, and could be used as tools to aid future theoretical work 
in resolving some of the complex questions in plasma physics. 

The results of Mahadevan (1998) encourages us to adopt the view that the protons 
in an ADAF, might in fact have a power-law distribution in energies. This gives rise to 
numerous testable predictions in the fluxes in various wavebands as well as the spectral 
shapes. First, all spectra from two body emission processes, in the optically thin limit, 
must have constant ratios in their fluxes. This is a result of all two body processes 
varying as the square of the density. The three processes in an ADAF, which depend 
on the square of the density, are bremsstrahlung emission and the creation of neutral 
and charged pions. From Figure ||b it is clear that these three processes always have 
the same relative amplitudes in the optically thin regime (v ^ 10 10 Hz). In particular 
the ratio of the gamma-ray to bremsstrahlung to e ± synchrotron fluxes are all the same, 
and are determined essentially by the ratios of their respective cross-sections and rates. 
While the cross-section for charged pion creation is greater than that for neutral pions, 
the synchrotron spectrum from the has a lower peak flux than the y-rays. This is 
due to the synchrotron radiation being emitted over a much broader frequency range. 

For low accretion rate systems, these emission processes dominate the radiation for 
frequencies v <: 10 12 Hz, which provides strong testable correlations among the fluxes in 
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the optical, X-rays, hard X-rays, and 7 -ray bands. Further, changes in the accretion 
rate, leads to changes in the the densities, which would result in simultaneous variations 
in the fluxes at these different frequencies. Observations of such variability would provide 
very interesting tests of the ADAF paradigm. It is apparent from Figure ^b that the 
best laboratories to verify these predictions would be in low luminosity AGN. 

I 11 addition to the predicted correlation in fluxes, there are also correlations in the 
expected spectral shapes and slopes. As discussed in §4, the ubiquitous presence of a 
hard power-law at ~MeV and ~ GeV energies provides two independent checks on the 
distribution of protons. I 11 an ADAF the detection of one requires the presence of the 
other. The power-law spectra presented here do not show any high energy cutoff. This 
is a direct consequence of postulating a power-law proton distribution which extends to 
very high energies. However, if we arbitrarily postulate a maximum proton energy (for 
example the energy at which the proton gyroradius equals the size of the ADAF) then 
the resulting spectra would also show a corresponding cutoff. 

The increased radio flux, from supermassive black holes, due to e ± radiation, facil¬ 
itates easier comparisons with observations, and can be used to constrain the sizes of 
ADAFs in AGN. Recently, Herrnstein et al. (1998) have compared observational upper 
limits to the 22 GHz emission from NGC 4258 with an ADAF model, and conclude 
that an ADAF cannot have radii larger than ~ 100 Schwarzschild. We might expect 
that the current results would produce more emission in the radio, and the size of the 
ADAF would therefore be further reduced. However, this is not the case since in order 
to produce more emission at 22 GHz, requires synchrotron radiation from e ± at radii 
r <: 100. I 11 fact the current results allows for a more definite upper limit to the size 
of the ADAF deduced by Herrnstein et al. (1998), since a larger outer radius implies 
more emission at 22 GHz from the e ± . This technique can also be used to constrain, or 
establish, the existence of ADAFs at the cores of nearby low luminosity ellipticals. 

Finally, we discuss two physical processes that have been neglected: (1) inverse 
Compton scattering of soft photons by the created e ± , and (2) the observational impor¬ 
tance of an annihilation line. First, the significance of inverse Compton scattering, by 
the created e ± , is determined by comparing the ratio of the total inverse Compton power 
to synchrotron power, Pcom/-Psyn = U^/Ub, where U p h, Ub are the photon and mag¬ 
netic field energy densities respectively. For an ADAF, U-^/Ub ~ ^emitted/-^advected 1, 
where T emitted , T advect ed are the emitted and advected luminosities, respectively, and 
inverse Compton scattering of the e ± is therefore not energetically important. 

Second, assessing the observational importance of an annihilation line requires a 
knowledge of the total energy in the line, as well as the line width. Since the annihilation 
cross-section is greatest for Lorentz factors 7 ~ 1, and most of the synchrotron emission 
occurs from electrons with 7 ~ 200 , the power in the annihilation line will be ~ P syn / 200 . 
In the case of Sgr A*this corresponds to ~ 10 33 erg s _1 . If the width of the line is narrow, 
then the line will be observable. However, the e ± that annihilate have thermalized 


with the ambient electrons at temperatures ~ 10 9 " 5 K. At such high temperatures, the 
annihilation line is sufficiently broadened (A E/E = \j2kT e /m e c 2 ~ 1 MeV), which 
results in a spectral line that is unobservable. 
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Fig. 1.— Figures showing the microphysics of e ± production, and the resulting syn¬ 
chrotron radiation. The left column is for a thermal distribution of protons where the 
curves are labeled by the dimensionless proton temperature, 6 P = kT p /m p c 2 . The right 
column is for a power-law distribution for three values of the energy index s. Figures 
(a) and (d) correspond the rate of creation of e ± for density of protons equal to unity. 
For a numberdensity n, the vertical axis must be multiplied by n 2 . Figures (b) and (e) 
show the resulting steady state distribution N( 7 ) of e ± . The y-axis plots C( 7 ) = 7 ^( 7 ) 
(cf. eq. 10). The synchrotron emission from the steady state distribution is shown in 
figures (c) and (f; cf. eq. 11 ). 
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Fig. 2.— Positron and electron synchrotron spectrum from an ADAF with m = 10 6 , 
m = 10~ 4 , for (a) a thermal distribution of protons, and (b) a power-law distribution 
with energy index s = 2.75. In the case of a power-law distribution, the spectrum which 
turns over at low (high) frequencies correspond to /3 a d v = 0.5 (0.9). 
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Fig. 3.— (a) Complete ADAF spectra (solid line) using a thermal distribution of protons 
and electrons, for different black hole masses and accretion rates. In all the calculations 
a = 0.3, /? a dv = 0.5. The e ± synchrotron spectrum (dotted line) is unobservable in all 
cases. For (m,rh) = (10,10“ 4 ), the 7 -ray spectrum and synchrotron spectrum have 
been multiplied by 10 ; the actual fluxes are ten times lower than shown. 
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Fig. 3.— (b) Complete ADAF spectra (solid lines) using a power-law distribution of 
protons and thermal electrons, for the same parameters given in figure 3a. Three spectra 
corresponding to proton energy indices s = 2.3, 2.75, and 3.3 are shown. The heavy 
solid line represents s = 2.75, with the spectra below and above corresponding to s = 3.3 
and s = 2.3 respectively. The dotted and dashed lines correspond to the emission from 
the e ± and the thermal electrons respectively. 
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Fig. 4.— The spectrum from an ADAF around Sgr A* for different values of the 
parameter A. The data are taken from Narayan et al. (1998b). A is the fraction 
of viscous energy that heats the protons into a power-law distribution. The solid, 
dotted, small dashed, dot-dashed, and long dashed lines (from high to low luminosity) 
corresponds to A = 1.0, 0.5, 0.25, 0.1, and, 0.0 respectively. The spectrum corresponding 
to A = 1 (solid line) is taken from Mahadevan (1998). 




